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Abstract. The existence of singularities of the solution for a class 
of Lax equations is investigated using a development of the fac- 
torization method first proposed by Semenov-Tian-Shansky and 
Reymann [TT], [5]. It is shown that the existence of a singularity 
at a point t = ti is directly related to the property that the ker- 
nel of a certain Toeplitz operator (whose symbol depends on t) be 
non-trivial. The investigation of this question involves the factor- 
ization on a Riemann surface of a scalar function closely related 
to the above-mentioned operator. An example is presented and 
the set of singularities is shown to coincide with the set obtained 
by classical methods. This comparison involves relating the two 
Riemann surfaces associated to the system by these methods. 

1. Introduction 

In this paper we investigate the existence of singularities of the so- 
lutions of Lax equations for a class of equations that applies to most 
finite-dimensional dynamical systems such as e.g. classical tops (see 
G.g. [3], [9], [10]). To that end we consider the time variable t to be 
a complex variable and determine the singularities of the solution in 
the complex plane. This is tied to the question of global existence of 
solutions for real t as the non-existence of singularities for real t implies 
global existence of the solution. Also, it is likely that full knowledge of 
the location of complex singularities may eventually give more insight 
into the dynamics of the system. 

Our approach is a development of the factorization method first pro- 
posed by Semenov-Tian-Shansky and Reymann [TT], [9] which in turn 
may be seen as a generalization of the AKS (Adler-Kostant-Symes) 
theorem that applies to finite dimensional algebras [TJ. 

To the best of our knowledge the first application of this method 
in the setting of an infinite-dimensional algebra appeared in [I], which 
focused on a restricted class of Lax equations. The absence in the liter- 
ature of more fully computed examples of application of this method is 
probably due to the fact that it involves the factorization of a continu- 
ous function on a contour in a Riemann surface (for a general treatment 
of this problem see [5]). 

Considering t as a complex variable and extending the class of Lax 
equations requires a new analysis of the results of [TJ, where we were 
able to avoid making some delicate assumptions like the differentiability 
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of the factors in the Wiener-Hopft factorization (see definition below) 
of the matrix function exp(£L ), where L Q is the value of the Lax matrix 
Lf at t = 0. This is done in Section|5]and continued in Section[3]for the 
question of location of the singularities of the solution. Our approach 
makes the treatment fully rigorous and, in our view, is crucial for the 
treatment in the context of a complex t. 

The main result is Theorem 13.1) which we state next. For this we 
note that the space [C^ {S 1 )]"" of Holder continuous nxn matrix-valued 
functions has a direct sum decomposition 

(i-i) pM(^)r=K] B ©[^]:. 

where [C^] n is the subspaceQof functions having an analytic extension 
to the unit disc D and [C~] is the subspace of functions admitting 
analytic extensions to C \ ID) that vanish at infinity. In what follows 

[qr] B 

Consider the Toeplitz operator 

(1.2) T G = P+GI + : [C+] B -> [C*]" 

where G = exp(£L ), 1+ is the identity operator on and P + 

is the projection of [C^] n onto [C^] associated to the decomposition 
(11. 2p . Theorem 13. II states that the Lax equation, 

(1-3) ^=[Lt + A ,L t ], 

has a solution in a neighbourhood of the point £j iff Tq is injective at the 
point ti. In the above equation A = PoL t where P : [C^S 1 )}™ n — > 
C nxn is a bounded linear operator. 

In the calculation of the singularities we need the notion of a Wiener- 
Hopf factorization. Let G : S 1 — > [C^\ nxn . We say that G possesses a 
Wiener-Hopf factorization (also called Riemann-Hilbert factorization 
and Birkhoff factorization [7], [S]) if G can be represented in the form 

(1.4) G = G-DG+, 



where G± and their inverses belong to [Cj] nX ™ and D = diag(r fcl , . . . , r k 
with ki > hi > ■ ■ ■ k n and r is a rational function with a zero in D and 
a pole in C \ ID). The factorization is said to be canonical if D = I n , 
where I n denotes the identity matrix. The above definition applies to 
functions belonging to other spaces (see e.g. [5]). In [C fl ] nxn G pos- 
sesses a Wiener-Hopf factorization (II. 3p iff G is invertible on S 1 . We 
recall from [6] that the operator Tq is invertible iff the factorization 
(11.4)1 is canonical. This is the basic result from operatior theory that 
will be used to locate the singularities of the solution of the Lax equa- 
tions. A direct consequence of Theorem 13.11 is Proposition 13.21 which 
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states that the solution of equation (11.31) has a singularity at t = t j iff 
the Riemann-Hilbert problem 

G$+ = $- (V g [c*] B , $- g [c^) 

has a nontrivial solution at i = tj. This is equivalent to saying that the 
Wiener- Hopf factorization of G is noncanonical for t = U. 

The paper ends with an example of a dynamical system that belongs 
to the standard Lax class considered in [I]. For this example it is 
possible to obtain the solution to the Lax equation by classical methods 
(integration of the system of ordinary differential equations) and thus 
obtaining its set of singularities. This enables us to compare it with 
the set of singularities derived by our method. A rather interesting 
point is that the classical approach and the Lax equation one lead to 
different Riemann surfaces. The two surfaces are closely related, as 
shown in Proposition 14.21 but the fact that they are different led us 
to derive several intermediate results in order to show that the sets of 
singularities obtained by the two approaches coincide. 

The study of the example given in Section 4 takes a large part of 
the paper but we believe that it not only illustrates the theory that 
we present here but also sheds some light into the relation between 
the classical methods and those based on the Lax equation - a point 
that may be obscure in the study of other finite-dimensional integrable 
systems, for example, in the study of some classical tops. 



2. Lax equation and Riemann-Hilbert problems 

In this section we generalize the results of |4, §2] replacing a neigh- 
borhood of the origin in the real variable t by a neighborhood of the 
origin in the now complex variable t (for what follows it is necessary 
to consider t as complex variable) and extending the class of equations 
considered. In [I] we studied a class of Lax equations of the form 

(2.D %=MM 

where the dynamical variables Lf , L t depend on a parameter A varying 
on the unit circle S 1 , L t is a matrix- valued Laurent polynomial in A 
and Lf is the part of L t analytic in the unit disc D. In [4] we called 
the above class the standard Lax class (it includes e.g. a special case 
of the Lagrange top). In this paper we study a class of Lax equations 
more general than the above one. It includes most finite-dimensional 
integrable systems. We write the equations of this class in the form 

(2-2) ^=[Lt + A ,L t ] 
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where Lf is defined as above and A = PoL t , with P being a bounded 
linear operator from the space [C 1 (f2)] nxn of matrix- valued Holder func- 
tions on S 1 to the space C nxn of constant matrix functions on S 1 (de- 
pending on t as a parameter). 

To state the first result in a rigorous way we need the definition that 
follows 

Definition 2.1. Let [C 1 (£l)] n><n be the space of continuously differen- 
tiable matrix functions with respect to t in a region flcC and define 
L t (X) G [C 1 (fi)] n n to be a Laurent polynomial of the form 

1 

(2.3) L t (X)= L? ] \ k (m G N, A G S 1 ), 

k=—m 

where l[® G C" xn . This gives for Lf(\) the expression 

(2.4) Lt(X) = J2Li k h k . 

Remark 2.2. In the case of the standard Lax class (Aq = in (12. 2p ) 
equation (12.21) together with formulas (12. 3p and (12. 4ft imply that 
is a constant of the dynamics. 

We can now state our first result which is a generalization of [H 
Theorem 2.3] extending the applicability of known formulas (see e.g. 
i) for h t . 

Theorem 2.3. Let L t be annxn matrix-valued function satisfying the 
Lax equation (12.21) in a simply- connected region Q containing the origin 
(in the variable t). Then L t is given in the region Q by the formulas 

(2.5) Lf = G j^LqG = G_^LqG^ 

where L = L t \ t=0 and G + , G_ satisfy in Q the linear differential equa- 
tions 

dC ~ dG ~ 

(2.6) -^± = (Lt + A )G + — = G^L--A ) 

subject to the initial conditions G + \t=o = G-\t=o = I n , where I n is the 
identity n x n matrix. 

Proof. The proof goes along the same lines as the proof of jH Theorem 
2.3] with Lf replaced by Lf ± A . We note only that the condition 
on the connectivity of fl is needed to ensure that equations (12.61) have 
well defined solutions throughout Q. □ 

Proposition 2.4. If the singularities of L t with respect to t are isolated 
then there exists a simply connected region Q = C \ B s where B s is the 
union of two cuts joining the singularities. 
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Proof. If the singularities are isolated we can denote them by t n with 
n G Z. Furthermore we can enumerate them in lexicographic order 
(t < t' iff Ret < Ret' or Ret = Ret' and Imt < hat'). Consider two 
consecutive points of this sequence, say tk,tk+i- Define one cut (B s +) 
as the union of line segments [t r _i, t r ] for r < k. Similarly, define the 
second cut (B s -) as the union of segments [i r ,t r+ i] for r > k. Then 
C \ (B s + U B s - ) is simply connected. □ 

Theorem 2.5. Let G = G-G+ where G-,G + satisfy equations (12. 6p 
in a simply connected region Q containing the origin (in the variable t) 
and the condition G\ t= o = I n . Then 

(2.7) G = exp(tLo). 

Proof. 

dG dG- ~ ~ dG+ 

G+ + Gr 



dt dt ' dt 
= G 



dt dt 



G A 



From equations (I2.6P 



Then ^ 

—— = G-L t G + = LqG, 
dt 

where we have introduced the expression for L resulting from the 
second of formulas (12. 5p . L = G-L t GZ l . Formula ( 12.71) now follows 
from the above equation. □ 

Theorem 2.6. 

(i) The factorization of G = exp(tLo), G = G-G+ is a canonical 
Wiener-Hopf factorization in the region Q of Theorem \2.5l 

(ii) Let G-G+ be another Wiener-Hopf factorization. Then 



(2.8) G + = G+F t , G_ = G-F t 



-i 
t ■ 



Proof. (i) Let G-G+ be the factorization of G obtained in Theo- 
rem [531 z.e. 

G = exp (tL ) =G-G+. 
This is a canonical Wiener-Hopf factorization of G in view of the 
properties of G-,G + resulting from equations (12.61) . 
(ii) Let be another (canonical) Wiener-Hopf factorization of 

G, obtained e.g. by solving a Riemann-Hilbert problem with co- 
efficient G ( C7$ + = $~). Then we have 

(2.9) G = G_G + = G_G + . 
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In the above relation both factorizations have factors that with 
their inverses are bounded analytic in their domains of existence, 
because it is assumed to be a Wiener-Hopf factorization 
and G + because the factors are assumed to satisfy equations 

From (12. 9p we have 

(2.10) GZ 1 G„ = G+Gf 

which implies that both sides equal a constant in A, but since we 
have a relation (I2.10p for each t, both sides must equal a function 
of t, independent of A. We write 

Kj Lx_ — 1_t_|_Lt_|_ — rf 

i.e. 

G— = G—Ff- , G-\. = FtG+. 

□ 

F t plays the role of a normalization factor at the point t for GL, G + . 
Theorem 2.7. Let F t satisfy the linear differential equation 

(2.11) § = A F t . 

Then the function L t given by 

(2.12) L t = F t L t F-\ 

with L t = G + L G^ 1 , satisfies the Lax equation (12. 2p . Here G + is as in 
Theorem \2.b\ (as are G_, G). 

Proof. We show that L t given by (I2.12p with F t satisfying (12. lip is a 
solution to equation (12. 2p . We have 

^ = AvF.LFr 1 + F.^-F,- 1 - F.L.F^AvF.Fr 1 

= lAoM + F t §F f K 

Since L t = G+LqG^ 1 it satisfies a standard Lax equation 

alL t 



(It 



L t , L t 



(see e.g. jU proof of Theorem 2.7]). Hence, noting that L t = F t L t F t , 
we get 

which is equation (12.21) . □ 
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Proposition 2.8. Equation (12. lip is equivalent to the equation in Q 

(2.13) f = FA 
where A = PoL t = PqG+LqG^ 1 . 

Proof. Let Aq = PoL t where Pq is a linear bounded operator as assumed 
in the definition of the right-hand side of equation (12.21) . Define 

(2.14) A = Ff l A F t . 

Noting that F t is independent of A, F t commutes with P Q which leads 
to 

Ao = PoFf L t F t = P L t = P G + L G + 1 . 
Introducing in equation (12.111) the definition (I2.14p of Aq we obtain 

^ = FtAoF^Ft = F t A , 
as required. □ 

Remark 2.9. Equation fl 2 . 1 3 [) is more convenient for calculating the 
solution of equation (12. 2p since Aq is known explicitly whereas Aq is 
not. 

Proposition 2.10. G = exp(tL ) has a canonical factorization at a 
point ti G f2 with factors differentiable w.r.t. t in a neighborhood of ti 
iff equation (12. 2p has a solution at the point ti. 

Proof. Sufficiency: Assume that equation (12.21) has a solution at a point 
ti G Q. Then by Theorem 12.31 there exist functions G_, G + satisfying 
(12. 6p in a neighborhood of ti (Q is open) which give the solution to 
equation (12. 2p . 

L t = G + LqG + 1 = GJ^LqG-. 

From Theorem 12.61 G + . G_ are related to the factors of another canon- 
ical factorization of G G + ) by the formulas 

G_ = G^Fi, Gr_|_ = Ft ^G+ 

where F t satisfies the differential equation (12. lip . Since C7_, G + and 
Ft are differentiable in a vicinity of ti it follows that the factors C7_, 
G + are differentiable too. 

Necessity: Assume that the factors C7_, G + of a canonical factorization 
of G are differentiable. Then 

dG r _dG dG + 



dt dt dt 

and letting L t = GZ 1 LqG- we get from the above relation 

dG+ n _ x 
dt 



L t '■- P+L t - —j— G + . 
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and, putting L t = F t L t F t 1 with F t satisfying (12. lip , we have (see the 
proof of Theorem 12.71) 

which is equation (12.21) . □ 



3. Singularities via the Riemann-Hilbert approach 

In this section we present our main result which enables us to locate 
the singularities of the solution to equation (12.21) without obtaining the 
explicit solution of the associated Riemann-Hilbert problem. Here the 
use of the factorization method is crucial since it allows us to translate 
the problem of the existence of singularities into an operator theory 
problem. 

We recall from the Introduction the direct sum decomposition of 

(3-i) K(&)] n =[ct] n ®[c;] n , 

where [C+] denotes the subspace of [C M (S 1 )} 71 of functions analytic in 
D and [C^*] is the subspace of analytic functions in C\D that vanish 
at infinity. We let P + : [C^(S 1 )] n — > denote the projection 

associated to this decomposition (so that kerP + = [C^]^)- 

Given a matrix G G [C M (S 1 )] nxn the corresponding multiplication 
operator in [C+] is denoted GI + . The composite 

is a Topelitz operator (with symbol G [HI Ch.l]) whose properties are 
closely related to those of its symbol G. We recall from [51 Ch.l] that 
the operator P + GI + is invertible iff G has a canonical Wiener-Hopf 
factorization 

G — Gr_(jr_j_, 

with(G ± ) ±1 e [CJOS 1 )]. 

We are now ready to state the main result of this section. 

Theorem 3.1. Let Tq be the Toeplitz operator P + GI + defined above, 
where G = exp(tL ). 

Then equation (12. 2 p has a solution in a neighborhood of a point t — ti 
iff the operator Tq is injective at that point, i.e., kerT^ is trivial. 

Proof. Sufficiency: We begin by proving that if kerT^ is trivial Tq 

is invertible. Firstly we note that G = exp(£L ) 6 C M (S 1 ) (in fact 

G G C°° {S 1 )) for every t G C Thus a factorization of G of the general 
form 

G = G-DG+, (Gi)* 1 G [Cj (^ x )] nx " , 
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where D is a diagonal nonsingular rational matrix, exists for all t G C 
(see e.g. Ch.l]). 

Since G G [C„ (S 1 )]™*™, det G G C M (S 1 ) and 

detG = exp((trL )t) ^ 0, 

for all A G S* 1 , it follows that Tq is Fredholm of zero index. This means 
that 

codimimTc = dim ker T^. 

It follows that, if ker Tq is trivial at t = ti, Tq is invertible at this point. 
It is easy to see that this is true in a neighborhood of ti. Hence G has 
a canonical factorization 

G = exp(tLo) = G-G+ 

in a neighborhood of t«. By Proposition 12.101 this implies that equa- 
tion (12 .2p has a solution at this point. 

Necessity: Assume that a solution to equation (I2.2p exists at t — t{. 
Then, by Proposition I2.10[ G possesses a canonical factorization at 
t ti , i. e. , 

G = G_C7 + , (G ± ) ±1 G [Cj (S 1 )]"^ . 

This is equivalent to the invertibility of Tq in a neighborhood of ti and 
thus kerTc is trivial. □ 

In the next two propositions we express the condition of Theorem 13. II 
in terms of the existence of solutions to a certain Riemann-Hilbert 
problem, which has the advantage of being easier to analyse. 



Proposition 3.2. Let To be the operator defined in Theorem \3. li Then 
ker To is nontrivial iff the Riemann-Hilbert problem 

with $~(oo) = 0, has non-trivial solutions. 

Proof, ker Tq being non-trivial means that the equation 

P+G$+ = 0, $+ G [C+] n 

has non-trivial solutions. Keeping in mind the direct sum decomposi- 
tion (I3.ip . we see that this is equivalent to saying that the Riemann- 



Hilbert problem in [C^ (5' 1 )] n 

G$ + = $- with $-(oo) = 0, 

has non-trivial solutions. □ 

Proposition 3.3. Letn = 2 in Propositio ns '.2 . Then the vector valued 
Riemann-Hilbert problem (on the Riemann sphere) 

(3.2) G$+ = $-, $-(oo) = 0, 
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given in Proposition 1 3. S\ is equivalent to a scalar Riemann-Hilbert prob- 
lem of the form 

(3.3) g^ + = 

on a compact Riemann surface E defined by the equation det(/ii2 — 
L(X)) = with subject to the condition 

(3.4) *~(oo x ) = 0, ^(oo 2 ) = 

where ooi, 002 are the poles of the meromorphic function given by the 
projection 

E -> P^C), (x,w) H> x 
(i.e., oo\, 002 are the points of E "at infinity"). 

Proof. It is proven in [5] that the Riemann-Hilbert problem (13. 2p is 
equivalent, for n = 2, to a scalar Riemann-Hilbert problem on E (13. 3p . 
The condition ( J3.4p is the translation of the condition $ _ (oo) = in 
( 13. 2 p to the Riemann surface. □ 

4. Example 

In this section we study a dynamical system for which the solution 
and, consequently, its singularities can be obtained by classical methods 
and compare the result obtained with that given by the method of 
Section [3J 

4.1. Dynamical system. We take the example presented in |4] which 
is given by the equations 



(4.1) 

where 

(4.2) 

with 



L t (X) 



v(X) 
w(X) 



u(X) 
-«(A) 



XeS 1 



(4.3) 



zX- 1 

aX + yX^ 1 + x, 

x 



v(X) 

u(X) = aX + yX^ 1 + x, a G C 
w(X) = aX + yX~ x 

and Lf being the polynomial part of L t (with respect to A). It can easily 
be seen that equation ( 14. ip together with ( 14. 2 p and (14. 3 p is equivalent 
to the following nonlinear system of differential equations 

dx dy dz 

4.4 — = -2az, -r = -2xz, — = 2xy 

y 1 dt dt dt y 

for the dynamical variables x, y, z. This system admits two integrals of 

the motion, namely, 



(4.5) 



A = x 2 - 2ay, 



B 



2 1 2 
y +z . 
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That these are invariants is easily checked by differentiating both sides 
of relations (14. 5 p and using equations f )4.4p . 

4.2. Classical solution. To obtain an equation of the movement in 
the variable x we begin with the first of equations (14.41) 

(xf = Aa 2 z 2 = Aa 2 (B - y 2 ) 

where x — Using relations (14. 5p yields 

( x f = Aa 2 B -{A- x 2 ) 2 = Aa 2 B - A 2 + 2Ax 2 - x 4 

from which we get 



(4.6) X = iy/f( 



X 



where p(x) = x A — 2Ax 2 + A 2 — 4a 2 B. The above equation means that 
(x, x) lies in an elliptic curve, i.e., the orbits of the dynamics lie in an 
elliptic Riemann surface. 

Before we integrate (14.61) we note that if we derive equations for the 
variables y, z we obtain equation (14.61) after an elementary transforma- 
tion on these variables as was to be expected. 

Integration of (14.61) gives 

r dx 

(4.7) / -4?= = it 



.to 



X 



where x is the value of x at t — 0, and the path of integration is 
understood to be on the Riemann surface £ defined by 

(4.8) w 2 = p(x) = (x 2 - x\)(x 2 - x 2 2 ) 
with the zeros of p(x), ±xi, ±X2, given by 

(4.9) x\ = A + 2aVB, x\ = A - 2aVB. 
It is useful to write (14.81) in the normalized form 

w 2 = xfxl(l — x 2 )(l — k 2 x 2 ) 

where x — xjx\ and k 2 = x\jx\. From now on we take as a definition 
of the Riemann surface £ the following equation 

(4.10) w 2 = (1 - x 2 ) (1 - k 2 x 2 ) , 

which corresponds to making the change of variables x h-> x/x±, w i— > 
w/(x\X2). With this notation, (14 .7p takes the form 



f jix _ , ^ |> 

Jx W{X) 
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with x = x/xi, x = xq/x\. It is convenient to write the above integral 
as a difference of two integrals as follows 

r dr r° dr 
(4.11) ix 2 t= to- *L (on U) 

Jo w(x) J w(x) 

We are looking for the singularities of the solution to equation ( 14. 6 p 
so we let x — > oo which leads to 

dots _ / f* _ / _ 

iK' or / — - iK' + 2K, 



(4.12) u 



lo w(x) Jo w(x) 

where K, K' are, respectively, the complete elliptic integral and the 
complementary complete elliptic integral of the first kind (see e.g. [2]). 
Hence from (14.111) . keeping in mind that (14. lip is an equation on E, 
we obtain 

ix 2 t = iK' -u + 4mK + 2iK', 

or 

ix 2 t = iK' + 2K - w + 4mK + 2iraK', n, m G Z, 

where 

^° dx 
w(x) 

The above formulas for t are equivalent to the single formula 

(4.13) ix 2 t = -uo + iK' + 2mK + 2mK', n, m e Z. 

This relation gives us the values of t at which singularities occur, i.e., 
where the solution blows up. 

4.3. Riemann-Hilbert solution. Next we derive a formula for the 
singularities of the solution to system (I4.4p using the method of Propo- 
sitions 13.21 and 13.31 To this end we need to formulate the Riemann- 
Hilbert problem (I3.2p for the function G = exp(tL ) in an associated 
Riemann surface. Recalling (13. 2p we have 

(4.14) G$ + = $~ with G = exp(tL ), 

and $ ± G [C^] 2 with the condition $~(oo) = 0. 
Taking into account that Lq can be diagonahzed as 

(4.15) L = SD S-\ 
with 



(4.16) S 



1 -1 

Lyi(A) <ji (A) 



(<?i(A) := a\ 2 + x X + yo) 



We identity E with the quotient of C by the lattice of periods of dx/ y/p{x). 
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and Dq = diag(A _1 /i, — A -1 ^) where /i = Xu with v given by the char- 
acteristic equation of L , 

(4.17) det (ul 2 - L (X)) = 0. 
From this equation we obtain 

(4.18) /i 2 = Pl (A) := a 2 A 4 - (x 2 - 2ay )X 2 + z 2 + y 2 
or, introducing the invariants A and B, 

(4.19) Pl {X) = a 2 X 4 -AX 2 + B. 

The explicit formulas for the zeros of Pi(A), ±Ai, ±A2, are given in 
(I4.33P below. Relation (I4.18P defines an elliptic Riemann surface, which 
is associated with L (or L t as it is independent of the dynamics). We 
denote by E\ the compact Riemann surface obtained by adding two 
points at infinity oo 1; oo 2 . 

Going back to (14.151) it follows from it that 

G = exp (tL ) = SDS- 1 

where D = diag(exp(£A _1 /i)> ex P( — ^"V))- 
Hence (I4.14p may be written as 

which, in terms of the components of $ ± , denoted (<pf, <pf), is written 
as 



(4.20) 



where d\ = exp(tA _1 /i); d 2 = exp(— tX~ l n). 

The above system is equivalent to the following single scalar equation 
(for more details see [3] or [S]) on a contour V that is the preimage of 
S 1 under the projection g: (A,/i) i— > A, 

(4-21) d(ri + Z -^ri)=^ + Z0 + ^- 



Qi J Qi 
Note that T has two connected components; we put d = d± on one of 
these components and al = d 2 on the other. In view of the expressions 
for d\ and d 2 we have 

(4.22) d = exp(^*), (A,/i)er. 

Concerning equation f)4.2ip . it is also useful to note that setting 
(72(A) = aA 2 — x X + y , we have 

(4.23) /i 2 -z 2 = gi (A)g2(A). 

It follows that, as a meromorphic function on Si, qi(X) has four zeros, 
two of which are zeros of \l + Zq and the other two are zeros of \i — Zq . 
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To solve (I4.2ip we note that d can be factorized on the Riemann 
surface as 

d = d_rd + , 

where (d + ) ±1 E (T) extends holomorphically to the preimage f2 + 
of D under the projection g and, similarly, (d - )^ 1 E C^(T) admits a 
holomorphic extension to the preimage Q~ of P(C) \ D. Finally, r is a 
rational function on E\. See [1] or [5] for more details. 

Note that all three factors in the above factorization depend on t. 
Introducing this factorization in (14.211) . we get 

(4.24) rd + ( 0+ + ^±^<fA = dZ 1 ( 2 - + ^-ft) = R, 



where R is a rational function on E\. 

For the computations that follow it is convenient to rewrite the Rie- 
mann surface E\ using the normalized equation 

(4.25) /i 2 = (l-A 2 )(l-^A 2 ), 

where k\ = X1/X2 and ±Ai,±A2 are the roots of pi(A) (cf. (I4.18P ) 
given in (14.331) below. This corresponds to making a change of variables 
A 1 — y A/Ai, \i 1 — y /x/(aAiA2). 

Also, from now on we identify E\ with its Jacobian, using the Abel 
map 

/■ (a ' m) dX 

(4.26) (X,»)^u= -, 

Jo V 

i.e., we consider all equations relating points of E\ as written on the 
quotient of C by the lattice of periods of the holomorphic form dX/n. 

Hence, keeping in mind that E\ is an elliptic Riemann surface (pi 
is a fourth degree polynomial), R can be expressed in elliptic theta 
functions. To this end we recall that we are solving (I4.24p with the 
conditions 0j r (oOj) = 0, for i,j = 1,2, where 001,002 are the two 
points at infinity □ in S 1 , which correspond to 00 under the projection 
g: (A,/i) I—)- A. Denoting by ip + ,ip~ the expression within parentheses 
in both sides of (14.241) . these conditions correspond to 

(4.27) V"(°oi) = 0, ^"(oo 2 ) = 0. 

Before introducing these conditions we note that, using the Jacobi theta 
function $1 that satisfies $1(0) = 0, R has the expression 

(4.28) *(„) = 7 'f-^f-^f-^. 

i)i{u - u )Vi{u - ui)Wi{u - u 2 ) 

where 7 E C and the zeros and poles of R are determined by the 
following conditions: 



We choose 00 1 such that /i <~ fciA 2 near 00 1. 
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(i) R has a pole at the point u$ corresponding to the pole of r in Q + 
(see |H Appendix B]); 

(ii) R has a zero at the point v o corresponding to the zero of r in + 
(see (U Appendix B]); 

(hi) R has two poles U\, u 2 at the zeros of q\ that do not coincide with 
the zeros of aX\X 2 \i + z (which is ji + z in (I4.24p written in the 
normalized coordinates of (I4.25P ); 

(iv) R has two zeros at points Vi, v 2 imposed by condition ( 14.27)) . i.e., 

V\ = OOi, v 2 = oo 2 . 

(v) The zeros and poles of R must satisfy Abel's condition: 

(4.29) vq — Uq = ui + u 2 — ooi — oo 2 (mod2iK' 1 + 4Ki), 

where Ki and K.[ are the complete elliptic and complementary 
elliptic integrals of the first kind of U\. 

From the analysis of the factorization of the function d given in (I4.22p 
(see [U Definition B.7 and Proposition B.9]) we obtain 

(4.30) v - u = 2at\ 2 

where X 2 is as in the expression for pi(A) (see text following (I4.19p ). 
Substitution of (I4.30P in (14.291) gives us the expression for the values of 
t for which singularities occur. Taking into account that ooi + 002 = 
2Ki + 2iK' x (mod 4mKi + 2mKi) we have 

(4.31) 2at\ 2 = Ul +u 2 + 2Ki + 4mKi + 2mK[ 

where U\,u 2 are the images under Abel's map of the zeros of q\ that 
do not coincide with zeros of a\i\ 2 fi + zq, i.e., 

f Xl/Xl dX f' X2/Xl dX 

ui= —, u 2 = —, 

Jo V Jo V 

where Ai,A 2 are de zeros of <?i(A) and Ai is a zero of Pi(A) given in 
(EOSD below. 

Thus (14.311) gives us the values of t leading to singularities of the 
solution of Lax equation (I2.2p as derived from the theory of Section [3j 

Remark 4.1. (i) Formula (I4.3ip was obtained without requiring an 
explicit formula for <fif, (fif corresponding to the factors of the 
canonical factorization of G, G = G_G + , although these functions 
can easily be obtained from (14.241) . replacing condition (iv) by the 
imposition of a zero at a chosen point v%. Then Abel's condition 
(v) gives the zero v 2 (see [1] for the details). The solution thus 
obtained gives the factors G_, G + of G providing t does not satisfy 
(I4.3ip . a result that could not be obtained in jl]. 
(ii) Formulas (I4.3ip and (I4.13P are not easily compared since they 
involve different Riemann surfaces. The appearance of distinct 
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surfaces when using different methods to study integrable systems 
is an intriguing phenomenon that occurs in other examples [3J. 

We show next that the two Riemann surfaces are closely related and 
that the two expressions for the singularities coincide. 

4.4. Comparison of solutions. We start by showing that the two 
Riemann surfaces S and E\ are related and this will enable us to 
express both formulas (14.13P and ( 14.31 j) on the same Riemann surface 
thus allowing for a comparison of the two results. 
The Riemann surface £ is defined by the equation 



,2 _ I -I ..-w'l 7.2, .2>i _ P( X 1 X ) 

\2 



W 2 = (l -x 2 )(l - fcV 



[XxX 2 ) 

withp(x) = x 4 — 2Ax 2 + A 2 — Aa 2 B = (x 2 — x\){x 2 — xf) and k — xi/x 2 , 
where 

(4.32) x\ = A + 2a^B, x\ = A - 2a\f~B. 

The Riemann surface E\ is defined by the equation 



,2 m \2 W -, , N /':(A|A) 

(aAiA 2 ^ 2 



/ x 2 = (l-A 2 )(l-A; 1 2 A 2 



^2) 

with Pl (A) = a 2 A 4 - AX 2 + B = a 2 (A 2 - A 2 ) (A 2 - A 2 ), where 

. A + y/A 2 - Aa 2 B - A - V 'A 2 - Aa 2 B 
(4.33) A? = - 2 , XI = - 2 . 

From the expression for p(x) and (I4.32p we have 

A 2 - Aa 2 B = xjxj 

Introducing these relations in (I4.33p gives 

,4.34) ^(^)\ A3=(-±- 

which leads to 

4.35 Ai = — , A 2 = — 

2a 2a 

where the sign in the square root is determined by a direct check on 

(I4.33p . From (14.351) we now get the relation between the moduli of the 

two surfaces 

(4.36) h = -r- = = 

A 2 x 2 + X\ 1 + k 

where k and k\ denote the elliptic moduli of E and Zi, respectively 
(see [2]). This shows that the surfaces are closely related as claimed at 
the end of Section 14.31 

Having obtained equality (I4.36P we are now in a position to state the 
following proposition relating S and Si. 
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Proposition 4.2. The following statements express the relation be- 
tween the Riemann surfaces E and E\: 

(i) For the elliptic moduli of E and E\, respectively k, k\, we have 

h = TTk 

(ii) There is a holomorphic map (p: E — > E\ given by 

E 3 (x,w) (A,/i) := (i(l + k)-, k X ~ l \ e E x . 

\ w w z J 

(Hi) Under the map (p of (ii) the points at infinity of E are mapped to 
Ox := (0,1) and (0, ±1), (±1,0), (±l//c, 0) are mapped to o 2 : = 
(0,-1), ooi and oo 2, respectively. 

(iv) The relation between the holomorphic forms of both surfaces is 
expressed by 

<P*[ — ) =-i(l + k — . 
V p. J w 

Proof, (i) was proven in (14.361) . 

The formula in (ii) is obtained by composing the two Gauss trans- 
formations corresponding in terms of elliptic moduli to k i— > k[ and 
k! i — y ki (see §39]). That it defines a map E — > E\ can be directly 
checked by a substitution of (I4.36P in /i 2 = (1 — A 2 )(l — k\}?). We note 
that this map is not injective; in fact it is 2 to 1. (Hi) is easily obtained 
by direct substitution in formula (ii). 

The expression (iv) follows directly from (ii) by differentiation. □ 



Before we attempt to formulate expression (I4.13P in E\ we are going 
to write u\ + u 2 of (14.311) as a single integral as in (I4.13P in order to 
make it possible to compare the two results. From (14.31 j) 



^ Jo V 
which we seek to write in the form 



/ ^ f dX 

(4.37) Ul + u 2 = —. 

Jo V 

To obtain £o we make use of the formula for the sum of arguments of 
the elliptic function sn (see [2]), 

lA ooX . . sn «i en «2 dn « 2 + sn « 2 cn Mi dn Mi 

(4.38) sn(tti + u 2 ) = -= — . 

1 — k\ sn z U\ sn z u 2 
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sn«! = — 



Ai 
V 



A 2 

snu 2 = — , 
Ai 



cn« 2 dnii 2 





1/2 






['-©'] 










-| 1/2 





Since q , i(AiA)g 2 (AiA) = (aAiA 2/ u) 2 — Zq (see ( I4.23P ) and Ai, A 2 are the 
zeros of q\ that are not zeros of aAiA 2 + z , we have 



(4.39) 
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where the factor l/(aAiA 2 ) comes from the normalization of fi. 
Using the above results in (I4.38p . taking into account, (14.371) gives 



(4.40) 



Ai + A 2 



Ai 1 - k\ \\ «AiA 2 

The denominator of the above formula can be simplified as follows 

\ 2 \ 2 \ 2 \ 2 
1 _ h 2 12 — i _ 12 

1 \4 ~~ 1 \2\2 

since A; 2 = Af/A|. Using f)4.39p we obtain 



X 2 X 2 
1_ Y| = i_ 
A X A 2 



a 2 A 2 A 2 B 



as A 2 A 2 = yo/a 2 , A 2 A 2 = 5/a 2 and B = y% + z 2 . 
Finally, 



(4.41) 



A^o qAiA 2 z 2 



5 _ x 
— — — A 2 

2 o 



where we have used the result Ai + A 2 = x /a. 

We shall now transform the terms on the right-hand side of (I4.13P 
into the surface S\. We first take the expression for uq given in (I4.12p 



(4.42) 

where 
(4.43) 



X °/ Xl dx 



W(X) 



X 



dX 

MA)' 



A = i(l + k) 



Xq/Xi 

w(x /xi)' 
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For the sake of simplicity in the calculations instead of trying to trans- 
form Ao into £o we prefer to take £o in ( I4.4ip and transform it as follows: 

i . . i\ j- . Xq lid 

4.44 £ = — A 2 = " A 2 

Zq w{Xq/Xi) XiX 2 

Xq 2ia X\ + X2 



w(xq/xi) X\X 2 2a 

= i^—^—Ai + k) = \Q. 

X 1 w{Xq/Xi) 

where we have used formula ( 14.451) and the relation w(xq/xi) = 2iazQ/(xiX 2 ) 
which is a consequence of (I4.4p and (14. 6 p for t = 0. (I4.35p . To end the 
calculation for the comparison of formulas (14.131) and (14.31 j) we need to 
derive relations between the complete elliptic integrals on S and Si. 
Using (ii) and (iv) of Proposition 14.21 we have 

r 1 dr 1 f°° d\ K' 

(4.45) K- ' dX - 1 dX - Kl 



o w(x) i(l + k) Jo n(X) i(l + k) 

(4.46) K'= T— = r — = 

1 1 Jq w(x) i{l + k)] 02 /i(A) i(l + *) 

Substitution of (Q4l) . (Q5]| and (051> in (j4TT31 now gives 

(4.47) - (1 + k)x 2 t = -i(l + A;)m + 2Ki + 4nK! + 2imK[. 

Using (I4.35P in (I4.3ip and the equalities k = xi/x 2 , \ 2 = (xi + x 2 )/2a, 
we see that formulas (I4.47P and (14.311) coincide. 
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